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Abstract 

The Belle experiment has used several methods to measure or constrain the an- 
gles (j)2, and ^3 (or /3, a, and 7) of the CKM unitarity triangle. The results 
are sm2(j)^ = 0.728 ±0.056 (stat) ±0.023 (syst) or (j)^ = (23.41^-^)° irom B° ^ 
J/il)K° decays (140 fb"!); (j)^ = (0-19)° or (71-180)° at 95.4% CL from 
7r+7r~ decays (253 fb^i); and ^3 = [68^j5(stat) ± 13 (syst) ±11 (model)]" 
from B^ ^ ,D°)K^ , ,D°) ^ tt+t:- decays (253 fb'^). These val- 
ues satisfy the triangle relation (f>i + (j)2 + (j)^ = 180° within their uncertainties. 
The angle is also determined from several b — + sqq penguin-dominated de- 
cay modes; the value obtained by taking a weighted average of the individual 
results differs from the i?° — > J/ip result by more than two standard devia- 
tions. The angle (/)2 is constrained by measuring a CP asymmetry in the decay 
time distribution; the asymmetry observed is large, and the difference in the 
yields of B^,B^ ^ tt+tt^ decays constitutes the first evidence for direct CP 
violation in the B system. 



1 Introduction 



The Standard Model predicts CP violation to occur in meson decays owing 
to a complex phase in the 3x3 Cabibbo-Kobayashi-Maskawa (CKM) mix- 
ing matrix This phase is illustrated by plotting the unitarity condition 
^ub^ud + ^cb^cd + ^tb^td = as vcctors in the complex plane: the phase re- 
sults in a triangle of nonzero height. Various measurements in the B system 
are sensitive to the internal angles (\)^, (^^i s-nd "^a (also known as /3, a, and 
7, respectively); these measurements allow us to determine the angles and 
check whether the triangle closes. Non-closure would indicate physics beyond 
the Standard Model. Here we present measurements of (/j^ and 02 obtained by 
measuring time-dependent CP asymmetries, and a measurement of (/)3 obtained 
by measuring an asymmetry in the Dalitz plot distribution of three-body de- 

cays. The results presented are from the Belle experiment 2) , which runs at the 
KEKB asymmetric-energy e+e~ collider operating at the T(4S') resonance. 

In Belle, pion and kaon tracks arc identified using information from time- 
of- flight counters, aerogel Cerenkov counters, and dE/dx information from the 

central tr acker ^) . B decays are identified using the "beam-constrained" mass 
= V^^boam - Pb and the energy difference AE = Eg- E^^^^^, where is 
the reconstructed B momentum, E^ is the reconstructed B energy, and E^^^^ 
is the beam energy, all evaluated in the e+e~ center-of-mass (CM) frame. A 

tagging algorithm ^) is used to identify the flavor at production of the de- 
caying B, i.e., whether it is B° or B°. This algorithm examines tracks not 
associated with the signal decay to identify the flavor of the non-signal B. The 
signal-side tracks arc fit for a decay vertex, and the tag-side tracks are fit for 
a separate decay vertex; the distance Az between vertices is to a very good 
approximation proportional to the time difference At between the B decays: 
Az ~ (d'ye) At, where /3j is the Lorcntz boost of the CM system. 

The dominant background is typically e^e^^qq continuum events, where 
q = u,d, s, c. In the CM frame such events tend to be jet-like, whereas BB 
events tend to be spherical. The sphericity of an event is usually quantified 

via Fox- Wolfram moments of the form = J2i j PiPj ^ei'^'^^^ij)' where i 
runs over all tracks on the tagging side and j runs over all tracks on either 
the tagging side or the signal side '''). The function Pg is the ith Legendre 
polynomial and 9^j is the angle between momenta p- and Pj in the CM frame. 
These moments are combined into a Fisher discriminant, and this is combined 
with the probability density hmction (PDF) for cosOq, where Og is the polar 
angle in the CM frame between the B direction and the z axis (nearly along the 
e~ beam direction). BB events are produced with a 1 — cos^O^ distribution 
while qq events are produced uniformly in cos^^. The PDFs for signal and qq 
background are obtained using MC simulation and M^^^-AE sidebands in data. 
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respectively. Wc use the products of the PDFs to calculate a signal likelihood 
£g and a continuum likelihood C^^ and require that Cg/{£,g + C^q) be above a 
threshold. 

The angles (f)^ and are determined by measuring the time dependence 
of decays to CP-eigenstates. This distribution is given by 

l-gAa; + g'(l-2a;)[^cos(AmAi)+5sin(ATOAt)] , (1) 

where q=+l (—1) corresponds to (-6°) tags, uj is the mistag probability, 
Aw is a possible difference in w between B^ and B^ tags, and Am is the B^-B^ 

mass difference. The CP-violating coefficients A and S are functions of the 
parameter A: ^ = (|A|2 - 1)/(|A|2 + 1) and S = 2Im(A)/(|A|2 + 1), where 

p A{B^^f) ^ VMi2 A{B^^f) \V,*^V,J A{B^^f) ' ^ ^ 

In this expression, q and p are the complex coefficients relating the flavor cigcn- 
states B^ and S ° to the mass eigenstates, M;^2 is the off-diagonal element of the 
B^-B ° mass matrix, and we assume that the off-diagonal element of the decay 
matrix is much smaller: ^ -^12- If only one weak phase enters the decay 
amplitude A(B° /), then \A(B° f)/A{B° ^ /)| = 1 and A = ry^e'^^, 
where rj^ = ±1 is the CP of the final state /. For the final states discussed 
here, |^| = or 



dN 
dAt 



(X e 



-At/r 



2 The angle 0^ 

This angle is most accurately measured using B^ J/ip K° decays^. The decay 
is dominated by a b^ccs tree amplitude and a b^scc penguin amplitude. The 

latter can be divided into two pieces: a piece with c and t in the loop that has 
the same weak phase as the tree amplitude, and a piece with u and t in the loop 
that has a different weak phase but is suppressed by sin^ 9q relative to the first 
piece. Due to this suppression, A{B'''^f) is governed by a single weak phase: 
Arg(V;^ V;* ). The ratio A(B° J/^j K^)/A{B^ J/tP K^) includes an extra 
factor {p/q)x = ^cd ^cs/ i^cd Ks) to account for the if ° oscillating to a Kg, and 

thus A = -[v,,v;j{v;,vj][v^,v,y{v:M][v:,vj{v^,v*^ = -e-'^1>K 

The CP asymmetry parameters are therefore S = sin 2(j)i, ^ = 0. To determine 
01, we fit the At distribution for S; the result is sin 2^^ = 0.728 ±0.056 (stat) ± 

iThis measurement includes B" J/ipK'^, J/ipKl, ^{2S)K'^, XcJ^s^ 
rj^K"g, and J/ijj K*° {K*° K°Tr°y, we use "S^ ^ J/ipK^" to denote aU 
six modes. 
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Table 1: Decay modes used to measure sin2(^i, the number of candidate events, 

the value of s\n2<f)i obtained, and the parameter A obtained [see Eq. (1)]. The 
J/tpK° result corresponds to 14O fb~^ of data; the other results corre- 
spond to 253 fb~^ of data. 



{CP) Mode 


Candidates 


sin2(/)]^ 


A 


{+)J/yjKl 


2285 
2332 


0.728 ± 0.056 ± 0.023 






139 ± 14 
36 ± 15 


0.08 ± 0.33 ± 0.09 


0.08 ± 0.22 ± 0.09 




398 ± 28 


0.74 ± 0.27 


-0.09 ±0.12 ±0.07 


(+ = 83%) 








(+)/o(980)Xg 


94 ± 14 


-OAT ± 0.41 ± 0.08 


-0.39 ±0.27 ±0.09 


i+)K'sK'sK's 


88 ± 13 


-1.26 ± 0.68 ±0.20 


0.54 ±0.34 ±0.09 




512 ±27 


0.65 ± 0.18 ± 0.04 


-0.19 ±0.11 ±0.05 




247 ± 25 


0.32 ± 0.61 ±0.13 


-0.11 ±0.20 ±0.09 




31 ±7 


0.76 ± 0.65 


0.27 ±0.48 ±0.15 



0.023 (syst), or (f)i = (23.4 ^2 4)° (^^ic smaller of the two solutions for (j>i). The 
fit result for A yields |A| = 1.007 ±0.041 (stat) ±0.033 (syst), in agreement with 

the theoretical expectation. These results correspond to 140 fb~^ of data 

There arc several decay modes that proceed exclusively via penguin am- 
plitudes (e.g., i? ° — > /sT " proceeding via b — » sss) or else are dominated by 
penguin amplitudes (e.g., — > {t]' /uj/w^)K^ proceeding via b — > sdd) but 
have the same weak phase as the b ccs tree amplitude. This is because the 
penguin loop factorizes into a c,t loop with the same weak phase and & u,t 
loop with a different weak phase; the latter, however, is suppressed by sin^ 9q 
relative to the former and plays a negligible role. We thus expect these decays 
to also have S = sm2(j)-^^, A = 0. There are small mode-dependent corrections 

(|A<S| < 0.10) to this prediction due to final-state rescattering Table 1 lists 

these modes and the corresponding values 

10) 

of sin 2(j)i obtained from fitting 
the At distributions; Fig. 1 shows these results in graphical form. Neglecting 
the small rescattering corrections and simply averaging the penguin- dominated 
values gives sin20j = 0.40 ± 0.13. This value differs from the S° J/tpK^ 
world average value by 2.4 standard deviations, which may be a statistical 
fluctuation or may indicate new physics. 
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Average 

b— >s penguin 



sin2(|)j(WA) 



"1 ' r 



I « I 



H 



H 

r 



0.08+0.33+0.09 

0.74+0.27to:i9 

-0.47±0.41±0.08 

0.65+0.18±0.04 

0.76+0.65to:l6 

0.32+0.61+0.13 

-1.26±0.68+0.20 

0.40+0.13 

0.726+0.037 



-1.5 



-0.9 -0.3 0.3 0.9 



1.5 



Figure 1: Values of sm2(f>i measured in decay modes dominated by b sqq 
penguin amplitudes, for 253 fb~^ of data. The average value differs from the 
world average (WA) value measured in B^^J/ijjK^ decays. 



3 The angle 4>^ 



This angle is measured by fitting the At distribution of ^tt+tt decays. The 
rate is dominated by a b^uud tree amphtude with a weak phase AigiV^^^ ^ud)- 
If only this phase were present, then A ^ [V,, V^JiV;, VJ] [V^, V:J{V:, V^,)] = 

duu penguin amplitude also 





e ■''2, and S — sin2(/)2, A = Q. However, a 6 
contributes, and, unlike the penguin in B'^ — > J/xpKg decays, the piece with 



a different weak phase is not CKM-suppressed relative to the piece with the 
same weak phase. The CP asymmetry parameters are therefore more compli- 
cated 11): 



1 

R 





p 








T 



2) sin S 



(3) 



c _ _ 





p 








T 



sm{(f>i — 4>2) cos S + sin 2(/)2 



sin 20]^ 



(4) 



5 



R=l-2 



COs((^l + 02 ) COS (5 



(5) 



where \P/T\ is the magnitude of the penguin amplitude relative to that of the 
tree amplitude, 5 is the strong phase difference between the two amplitudes, 

and is known from J /ipK^ decays. Since Eqs. (3) and (4) have three 

unknown parameters, measuring A^^ and S^^ determines a volume m. 5'\P/T\- 
<p2 space. 

The most recent Belle measurement 

12) 

uses 253 fb ^ of data; the event 
sample consists of 666 ± 43 5° — > tt+tt" candidates after background sub- 
traction. These events are subjected to an unbinned maximum likelihood 
(ML) fit for At; the results are A^^ = 0.56 ± 0.12 (stat) ± 0.06 (syst) and 
<S^^ = — 0.67 ± 0.16 (stat) ± 0.06 (syst), which together indicate large CP vio- 
lation. The nonzero value for A^^^ indicates direct CP violation. Fig. 2 shows 
the At distributions for q = ±1 tagged events along with projections of the ML 
fit; a clear difference is seen between the fit results. 

The values of and S^^ determine a 95.4% CL (2(t) volume in d-\P/T\- 
<p2 space. Projecting this volume onto the S-\P/T\ axes gives the region shown 
in Fig. 3; from this region we obtain the constraints |P/T| >0.17 for any value 
of 6, and -180° <(5< -4° for any value of \P/T\. 

The dependence upon 5 and \P/T\ can be removed by performing an 

isospin analysis 

13) 

of _B— s-TTTT decays. This method uses the measured branch- 
ing fractions for B tt^tt^ , tt^tt", tt^tt" and the CP asymmetry parameters 
.4 + S . _ , and ^ o o • We scan values of 6^ from 0° — 180° and for each 
value construct a % based on the difference between the predicted values for 
the six observables and the measured values. We convert this iiito a con- 
fidence level (CL) by subtracting off the minimum value and inserting the 
result into the cumulative distribution function for the distribution for one 
degree of freedom. The resulting fimction 1 — CL is plotted in Fig. 4. From 
this plot we read off a 95.4% CL interval 02 = (0-19)° or (71-180)°, i.e., we 
exclude the range 20° -70°. 



4 The angle (p^ 

The angle (p^ is challenging to measure by fitting the At distribution, as the two 
requisite interfering amplitudes have very different magnitudes, and the small 

ratio of magnitudes multiplies the 03-dependent term [sin(20]^+03± S)] 14). As 
an alternative, one can probe cp^ via interference in the Dalitz plot distribution 
of — !■ {D^ , D^)K^ decays: the additional phase (f>^ causes a difference 

between the interference pattern for B^ decays and that for B~ decays 

15). 
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Figure 2: The At distribution of background- subtracted B'^ , B'^ —ht^tt candi- 
dates (top), and the resulting CP asymmetry [N(B°)-N{B°)]/[N(B°y\-N{B°)] 
(bottom). The smooth curves are projections of the unbinned ML fit. 

We study this asymmetry by reconstructing B^ {D^ ,D^)K^ decays 
in which the or decays to the common final state _ftr° tt+tt". Denoting 
m{Kg,T:+) = TO+, miKg,-!:') = m_, ivT^ tt+tt^) = A{m+,m~), and 

^ TT+TT^) = A(m^ ,m^) = A(m^,r7i+) (i.e., assuming CP conserva- 
tion in decays), we have 

(6) 

KItt+t:-) = A{mt,mX) + re''^^-'^^^A(mX,m\), 

(7) 

where denotes + D"), r is the ratio of magnitudes of the two ampH- 
tudes \A{B+ ^D^K+)/A{B+ ^T)^K+)\, and 5 is the strong phase difference 



A{B+^D°K+, & 
A{B- ^D°K-, & 
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Figure 3: Projection of the 68.3% CL (dashed) and 95.4% CL (solid) volumes 
in 5-\P /T\-(j)2 space onto the 5-\P/T\ axes. From the solid contour we obtain 
the constraints |P/T|>0.17 and -180°<,5<-4° (95.4% CL). 

between the amplitudes. The decay rates are given by 



A{B+ ^D°K+-^ iK°s TT+TT- ) A'^ 



\A{m\,mi)\'^ +r'^\A{m'i,m.\)\'^ + 



2r\A{ml,ml)\\A(m'i,ml)\cos{S + (1)3 + 9) (8) 



A{B- D°K- {K°s TT+n-JK- 



r''\A{ml,m'i)f + \A{m'i,ml)\'' + 



2r\A{ml,m'i)\\A{nvi,ml)\ cos{S - ^3 + 6») , (9) 
where 6 is the phase difference between A{m\ , m?_) and A(w?_ , m\_) and varies 



over the Dalitz plot. Thus, given a 



decay model A(m 



one can fit the Dalitz plots to Eqs. (8) and (9) to determine the parameters 
r, 6, and (p^. The decay model is determined from data, i.e., D*^ — > i^^ tt+tt" 
decays produced via e~^e~^cc. 

The data sample used consists of 253 fb~^; there are 209 ±16 B^^D^K^ 
candidates with 75% purity, and an additional 58 ± 8 B^^D°*K^ {D°* 

D^tt'^) candidates with 87% purity , The background is dominated by qq 
continuum events in which a real is combined with a random kaon, and 
random combinations of tracks in continuum events. The Dalitz plots for the 
final samples are shown in Fig. 5. 
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Figure 4: The result of fitting the branching fractions for B -k^i:~ , tt^tt^, 
TT^TT^ and the CP asymmetry parameters A^+^- , iS^+^_, andA^o^o, as a func- 
tion 0/^2 T'^ee text). The vertical axis is one minus the confidence level. The 
horizontal line at 1 — CL = 0.046 corresponds to a 95.4% CL interval for 02- 

The events are subjected to an unbinned ML fit for r, 6, and The 
decay model is a coherent sum of two-body ampHtudes and a constant term 
for the nonresonant contribution: 

N 

A{ml , ) = ^ a^. e'"^ Aj {ml ,m^_) + a„„„,,, e'"— , (10) 

where a^, a^, and Aj are the magnitude, phase, and matrix element, respec- 
tively, of resonance j; and A^= 18 resonances are considered. The parameters 
and aj are determined by fitting a large sample of continuum Kg tt+tt" 

decays. The dominant intermediate modes as determined from the frac- 
tion / \a^Aj P dm,l dm^_/ J \A{m\ ,m'i)\'^ dm\ dml are ii'*(892)+7r- (61.2%), 
i^gpO (21.6%), nonresonant X^TT+Tr" (9.7%), and i4:^(1430)+7r- (7.4%). 

The central values obtained by the fit are r = 0.25, 6= 157°, and c/ig =64° 
for B+-^D°K+; and r = 0.25, (5 = 321°, and cf>^ = 75° for B+^D*°K+. The 
errors obtained by the fit correspond to Gaussian-shaped likelihood distribu- 
tions, and for this analysis the distributions are non-Gaussian. We therefore 
use a frequentist MG method to evaluate the statistical errors. We first obtain 
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Table 2: Results of the Dalitz plot analysis for r, S, and (f>^. The first error 

listed is statistical and is obtained from a frequentist MC method (see text); 
the second error listed is systematic but does not include uncertainty from the 
-D°— s-JCgTr+TT" decay model; the third error listed is due to the decay model. 



Parameter 


B+^D°K+ B+^D*°K+ 


r 
S 

03 


0.21 ± 0.08 ± 0.03 ± 0.04 0.12 1° ± 0.02 ± 0.04 
157° ± 19° ± 11° ± 21° 321° ± 57° ± 11° ± 21° 
64° ± 19° ± 13° ± 11° 75° ± 57° ± 11° ± 11° 



a PDF for the fitted parameters r, 5, ^3 as a function of the true parameters 
f, S, ^3. We do this by generating several hundred experiments for a given set 
of f, S, 4>^ values, with each experiment having the same number of events as 
the data, and fitting these experiments as done for the data. The resulting 
distributions for a_|_ = r cos((5 ± ^3) and f3±=r sm{S ± ^3) are modeled as Gaus- 
sians (G) with mean values a^. and /3_|_ and common standard deviation a, and 
the product G{a_f_ — a^) ■ G{a_— a_) ■ G(/3^— /3^) ■ G{(3_ — (3_) is used to ob- 
tain the PDF 'P{r, 5, (jy^lf, S, (j)^). With this PDF we calculate the confidence 
level for {f, 5, ^3} given the fit values {0.25, 157°, 64°} for B+^b^K+ and 
{0.25, 321°, 75°} for ^b*^K+ . The resulting confidence regions for pairs 
of parameters are shown in Fig. 6. The plots show 20%, 74%, and 97% CL 
regions, which correspond to one, two, and three standard deviations, respec- 
tively, for a three-dimensional Gaussian distribution. The 20% CL regions are 
taken as the statistical errors; the values that maximize the PDF are taken as 
the central values. Of the two possible solutions {5, ^3) or [5 -|- tt, (^3 -|- tt), we 
choose the one that satisfies 0° < ^3 < 180°. 

All results are listed in Table 2. The second error listed is systematic 
and results mostly from uncertainty in the background Dalitz plot density, 
variations in efficiency, the m^^ resolution, and possible fitting bias. The third 
error listed results from uncertainty in the — >ii'^7r+7r~ decay model, e.g., 
from the choice of form factors used for the intermediate resonances and the 
dependence of the resonance widths. 

We combine the B+ ^b^K+ and B+ ^b*°K+ results by multiplying 
together their respective PDF's, taking the parameter (f>^ to be common be- 
tween them. This gives a PDF for the six measured parameters r^^, 5^, 4'z(i)^ 
Tj, ^2> '^3(2) terms of the five true parameters f^, 5^, fj, ^3. The value 
of (^3 that maximizes the PDF is taken as the central value, and the 3.7% 
CL interval prescribed by the PDF (corresponding to Icr for a five- dimensional 
Gaussian distribution) is taken as the statistical error. The systematic error 
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is taken from the 5+ D^K'^ nieasurcmcnt, as this sample dominates the 
combined measurement. The overall result is 

cj)^ = [68 t\t (stat) ± 13 (syst) ±11 (decay model) ] ° . 

The 2cr confidence interval including the systematic error and decay model 
error is 22° < (^ig < 113°. 

In summary, the Belle experiment has measured or constrained the an- 
gles 4>2, and 03 of the CKM unitarity triangle. We obtain sin 20 = 
0.728 ± 0.056 (stat) ± 0.023 (syst) or 0i = (23.4 12-7^ ° ^ith 140 ft-^ of data; 
02 = (0-19)° or (71-180)° at 95.4% CL with 253 fb"! of data; and 03 = 
[ 68 t\t (stat) ± 13 (syst) ±11 (decay model) ] ° with 253 fb"! of data. Within 
their uncertainties, these values satisfy the triangle relation 0]^ +02 + 03 ~ 180°. 
The angle (j)^ is measured from i?" — > J /%p decays and also from sev- 
eral b sqq penguin-dominated decay modes; the value obtained from the 
penguin modes differs from the — » J/tpK^ result by 2.4 cr. The 02 con- 
straint results from measuring the CP asymmetry coefficients A^^^ and S^^ 
in 5° ^TT+TT" decays; the results are =0.56 ± 0.12 (stat) ± 0.06 (syst) 
and S^^ = —0.67 ± 0.16 (stat) ± 0.06 (syst), which together indicate large CP 
violation. The nonzero value for A^^^ indicates direct CP violation; the statis- 
tical significance (including systematic uncertainty) is 4.0 a. These values also 
imply that the magnitude of the penguin amplitude relative to that of the tree 
amplitude {\P/T\) is greater than 0.17 at 95.4% CL, and that the strong phase 
difference {S) lies in the range (—180°, —4°) at 95.4% CL. The 03 measurement 
is obtained from a Dalitz plot analysis of B^^D^*^^K^, D^-^K^n+Tr- de- 
cays; the statistical significance of the observed (direct) CP violation is 98%. 
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Figure 5: Dalitz plots of JC^tt+tt decays obtained from samples of 

B+ D°K+ (top left), B- D°K- (top right), B+ D*°K+ (bottom 
left), and B- ^D*°K- (bottom right). 
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Figure 6: Confidence regions for pairs of parameters: the left-most plots cor- 
respond to r-(p^ and the right-most plots to 6-(j)^. The top row corresponds to 
B^^D^K^ decays and the bottom row to B^^D*°K^ decays. 
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